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(3 , ABSTRACT 

■ We present a method to simulate the polarization properties of extragalactic radio 

I sources at microwave frequencies. Polarization measurements of nearly 2 x 10^ sources 

CO . at 1.4 GHz are provided by the NVSS survey. Using this catalogue and the GB6 survey, 

' we study the distribution of the polarization degree of both steep- and flat-spectrum 

fH I sources. We find that the polarization degree is anti-correlated with the flux density 

for the former population, while no correlation is detected for the latter. The available 
high-frequency data are exploited to determine the frequency dependence of the distri- 
bution of polarization degrees. Using such information and the evolutionary model by 
IXoffolatti et alJ ((1.998), we estimate the polarization power spectrum of extragalactic 
radio sources at ^ 30 GHz and their contamination of CMB polarization maps. Two 
distinct methods to compute point-source polarization spectra are presented, extend- 
ing and improving the one generally used in previous analyses. While extragalactic 
radio sources can significantly contaminate the CMB E-mode power spectrum only 
at low frequencies {v < 30 GHz) , they can severely constrain the detectability of the 



o 



■ CMB B-mode up to ~ 100 GHz. 
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1 INTRODUCTION 

In the last year, two experiments provided the first ob- 
servational evidence of the polarization of the Cosmic Mi- 
crowa ve Background (CMB) radiation: DASI iKovac et alJ 
I 2OO3) achieved a direct me asu re of E-mode polarization ( see 
iKamionkowski et al .119971 and lZaldarriaga fc Seliakll997lfor 
a defin ition of E and B modes) . while WMAP JseimB^et alJ 
detected the cross-correlation between the CMB 
temperature anisotropies and E-mode polariz ation. Riving 
an es timate of the reionization optical depth llKoeut et alJ 
l2003l) . A large number of experiment s have been planned 
to measure the CMB polarization (seelCecch ini et alll2002t 
FHananv fc Olive 2003), motivated by its huge information 
content. Besides probing the ionization history, the CMB po- 
larization can open a window on the primordial phase of the 
Universe through the detection of the B-mode, i.e. the curl 
component of the polarization field. In fact, the B-mode po- 



larization is generated by tensor metric perturbations and, 
in inflation models, its amplitude is directly proportional to 
the energy scale at which the inflation has occurred. How- 
ever, the detection of this component is a real challenge, not 
only due to the very low level of the signal, but also to the 
presence of foregrounds and effects that can mix the E- and 
B-modes of the CMB polarization. For example, the gravi- 
tational lensing produced by large scale structures converts 
a fraction of the CMB E-mode component to the B-mode 
one. Although the effect is of a few percent or below, the 
lensi ng-induced curl mode can overwhelm the pri mordial 
one dZaldarriaea fc SebakllTool iKnox fc Sonell2002^ . 

Our knowledge of foreground polarization is still poor. 
Among the Galactic emissions which dominate foreground 
intensity fluctuations on large angular scales, free-free is 
nearly unpolarized while synchrotron is highly polarized; 
its polari zation properties have b een studied at GHz fre- 
quencies llTucci et alJ I200I 2002; iBacciealupi et alJ 1200 ll: 
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iBruscoli et al .120021) ■ but we have little or no direct informa- 
tion at frequencies of interest for CMB studies. At mm wave- 
lengths, the Galactic polarized signal is due essentially to the 
dust emission. Estimates of this contribution have been so 
far pr ovided only by models based on H I maps jPrunet et alJ 
Il998r or on the starlight polarization llFosalba et al.ll2002D . 
although direct sub-mm measurements over a substantial 
fraction of the sky have been recently pr ovided by the 
ARCHEOPS experiment (iBenoit et alji2003i) . 

In this paper we deal with extragalactic radio sources, 
which are expected to be the dominant polarized foreground 
on small angular scales. Accurate studies of their contribu- 
tion to intensity fluctuations, based on evolutionary models 
fitting the available data at many frequen cies, have been 
publis hed. In our analysis we adopt the iToffolatti et alJ 
il99g h (hereafter T98) model in order to predict the num- 
ber counts of extragalactic radio sources at cm and mm 
wavelengths. This model reproduces very well the number 
counts of different classes of bright radio sources at GHz fre- 
quencies (see Sections 121 and I^J ; moreover, the accuracy of 
its high-frequencies predictions has been confirmed by the 
recent VSA (Taylo r et al. 2001; Waldram et al 200^), CBI 
jMason et alj|2003l) and WMA P jBennett et alj|2003bl) sur- 
veys between 15 and 44 GHz jToffolatti et alJl2003ll More 
specifically, the WMAP survey has detected 208 sources in a 
sky area of 10.38 sr at |6-^^| > 10° at flux densities S > 0.9- 
1.0 Jy. The T98 model predicts 270-280 sources in the same 
sky area, so that the average offset is a factor of ^ 0.75. 
Moreover the distribution of spectral indexes in the WMAP 
sample peaks around a = 0.0, which is exactly the mean 
spectral index for fiat-spectrum sources adopted by T98, 
and the fraction of steep-spectrum (i.e., a > 0.5, S oc 
sources is of ^ 12%, to be compared with a predicted frac- 
tion of ~ 10%. As a curiosity, we also note that the brightest 
source detected by WMAP has a fiux density of 5' ~ 25 Jy 
which corresponds exactly to the value for which the model 
predicts 1 source all over the sky. 

Most of the available information on polarization prop- 
erties of extragalactic sources refers to GHz frequencies, 
while for CMB studies we are interested in frequencies > 
30 GHz. There are at least two reasons to expect polarization 
properties to be frequency dependent, and, in particular, 
higher polarization at higher frequencies. First, at low fre- 
quencies substantial depolarization may be induced by Fara- 
day rotation. Determinations of rotation measures for radio 
galaxies and quasars yield values from tens to lO^radm"^ 
JO'Deall989l : lTavloil2000l:lPentericci et alj200ol : lMesa et all 
|20oJlienceforth M02), implying that depolarization may be 
significant up to ~ 10 GHz. Second, especially in compact 
objects, as the observing frequencies increase emitting re- 
gions are closer and closer to the nucleus, where the degree 
of the ordering of magnetic fields and, as a consequence, the 
polarization degree, may be higher and higher. 

Preliminary estimates of the contribution of extragalac- 
tic radio sources to the pol arization angula r powe r spectrum 
(APS) were worked o ut bv lDe Zotti et al.l ^1999^ and M02. 
iDe Zotti et al.1 il999l) simply scaled the temperature power 
spectrum by the mean squared polarization degree of ra- 
dio sources, assuming a Poisson spatial distribution. M02, 
instead, co mputed the APS dire ctly from the NVSS polari- 
metric data jGondon et alll998D . estimating a factor of 3 in- 
crease in the polarization degree of steep-spectrum sources 
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Figure 1. Differential number counts, normalized to S^/^, of 
steep- (upper data points) and flat-spectrum (lower data points) 
NVSS sources, compared with predictions of the T98 model 
(steep-spectrum sources: short dashed curve; flat-spectrum 
sources: dotted curve; "starburst" galaxies: long dashed curve). 
The upper part of error bars for steep spectrum sources allows 
for the possibility that NVSS sources missing a GB6 counterpart 
belong to this population. The dots show the total differential 
counts in terms of polarized flux. 

from 1.4 GHz to ^ 20 GHz and assuming that, for fiat- 
spectrum sources, the polarization degree either is frequency 
independent or increases by a factor of 3. 

In this paper, we improve on previous estimates by 
studying the probability distribution of the polarization 
degree for both steep- and flat-spectrum sources at i/ ^ 
15 GHz, providing a recipe to extrapolate the distributions 
to higher frequencies. Our results consist of a technique to 
generate simulated catalogues of extragalactic radio sources 
and also in accurate expressions to estimate the polarization 
APS, after having removed the strongest polarized sources. 



2 POLARIZATION PROPERTIES OF NVSS 
SOURCES 

The starting p oint of the analysis is the NRAO VLA Sky 
Survey (NVSS JCondon et alJll998l) . This survey covers Q. ~ 
10.3 sr of the sky with S ^ -40° at 1.4 GHz, providing 
the fiux density 5* and the Stokes parameters Q and U of 
almost 2x 10^ discrete sources with S > 2.5 mjy. The images 
obtained from interferometric data have 8 = 45" FWHM 
resolution and the rms brightness fiuctuation in the Q and 
U parameters is cr ~ 0.29 mJy beam" ^. 

Information on the spectral index a of NVSS sources is 
obtained exploitin g the Green Bank 4.85 GHz survey (GB6, 
iGreeorv et alll 996'). which covers the declinations 0° < 5 < 
75° (O ~ 6.07 sr) to a flux limit of 54.85 = 18mJy with a 
resolution FWHM= 3.5 arcmin. We have cross-matched the 
positions of NVSS sources with Si. 4 ^ 100 mJy with the 
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Table 1. Polarization degree at 1.4 GHz of steep— and fiat— spectrum sources. 



flux(mjy) 


N 


N(T98) 


fiat 
Nn<i%(%) 


median(%) 


mean(%) 


N 


N(T98) 


steep 
Nn<i%(%) 


median(%) 


mean(%) 


100-200 


2305 


2144 


43.3 


1.33 


2.16 


15621 


15409 


35.8 


1.77 


2.70 


200-400 


980 


1045 


43.5 


1.24 


2.01 


6925 


6006 


40.0 


1.52 


2.40 


400-800 


473 


486 


40.2 


1.50 


1.95 


2552 


2428 


41.2 


1.44 


2.25 


> 800 


260 


364 


45.7 


1.12 


1.92 


1096 


1165 


47.0 


1.14 


2.02 



GB6 catalogue, taking all matches with position oflFsets less 
than 3 times the uncertainty on the GB6 position. Sources 
with Galactic latitude |6| ^ 2° have been rejected, in order to 
guarantee that nearly all the sources are extragalactic. The 
100 mjy flux limit has been chosen to have > 3cr detections 
of polarization down to a percent level. About 86% of NVSS 
sources turned out to have a GB6 counterpart. Whenever 
more than one NVSS source falls within the GB6 beam, as 
a consequence of the better NVSS angular resolution which 
may lead to individually resolve multiple components, we 
have summed up their fluxes, corrected for the effect of the 
GB6 beam. We ended up with a sub-sample of 29299 NVSS 
sources, of which ~ 87% are steep-spectrum [a > 0.5) and 
~ 13% are flat-spectrum (a ^ 0.5). 

In Fig. we show the differential number counts of 
NVSS sources a s a functi on of the total and linearly po- 
larized {PI — \/QM-T7^) flux density. The predictions of 
the T98 model (for a flat CDM cosmology with = 0.7 
and h — 0.7) are also displayed. According to the model, the 
dominant contribution at high and moderate fluxes comes 
from steep-spectrum sources, while "starburst" galaxies are 
negligible for S above a few tens of mJy. The fraction of 
flat -spectrum sources increases with flux density and is im- 
portant at the Jy level. As shown by Table and by Fig. 
Q the observed number of flat- and steep-spectrum sources 
in total flux-density bins compares favourably with the pre- 
dictions, N(T98), of the T98 model. The latter only weakly 
overestimates the observed number of flat-spectrum sources 
for S > 800 mJy; such discrepancy will be taken into account 
when computing their contribution to the polarization power 
spectrum in Section |3 

It is interesting to note that the normalized PI number 
counts keep nearly constant over all the flux density range 
where they are defined. The different shape compared with 
the total power number counts indicates a dependence of 
the polarization degree on flux density (see below) . 

Following the analysis of M02, we study the distribution 
of the percentage polarization degree (11 = 100P//5') at 1.4 
GHz for steep- and flat-spectrum sources for different flux- 
density intervals (see Figs. |5| and |^. Table reports the 
median and the mean value of 11 for such distributions and 
the percentage of sources with 11 lower than 1% (Nn<i%). 
The low-II tail of distributions is contaminated by noise 
and residual instrumental polarization; however, while the 
relevance of the noise decreases when the polarized intensity 
increase s, the latte r is proportional to the flux density of 
sources. lGondon"et al. ( 1998) found that the instrumental 
polarization is ~ 0.12% of S for a large sample of sources 
stronger than 1 Jy and estimated that, in any case, it should 
be less than 0.3%. 

The results of Tabled highlight an anti-correlation be- 
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Figure 2. Polarization degree distribution of steep-spectrum 
NVSS sources for the flux density intervals specified in each panel. 
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Figure 3. As in Figure|2]but for flat-spectrum sources. 
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Table 2. Results of the test for 11 distributions at dif- 
ferent fluxes. The distributions are binned in intervals of the 
1% width. Probability values above the diagonal refer to flat- 
spectrum sources, those below the diagonal, to steep-spectrum 
sources. 



(mjy) 


100-200 


200-400 


400-800 


> 800 


100-200 




0.28 


0.91 


0.95 


200-400 


2.1 X 10-4 




0.12 


0.06 


400-800 


1.8 X 10"'' 


4.8 X 10-2 




0.52 


> 800 


1.4 X lO"'' 


6.1 X 10-3 


0.16 





tween the polarization degree and the flux density for steep- 
spectrum sources: the median value of D steadily decreases 
from 1.8% at 100-200 mJy to 1.1% at S ^ 800 mJy (we con- 
sider the median because it is less affected than the mean by 
instrumental polarization). A similar trend is not found for 
flat-spectrum sources, whose median values show only small 
fluctuations, all compatible with a constant value at the 1-cr 
level. Pearson's linear correlation coefficient r between the 
flux density and the polarization degree of sources is low 
for both flat (-0.017) and steep-spectrum (-0.034) sources, 
yielding a 30% probability of the null hypothesis (i.e., no cor- 
relation) in the flat case. However, the test clearly rules 
out the possibility that the polarization degree distributions 
of steep-spectrum sources in the range 100-200 mJy and at 
higher flux density are drawn from the same parent distribu- 
tion (probabilities of the order or less than 10"''; see Table 
On the contrary, in the case of flat -spectrum sources, the 
test indicates consistency with the same parent distribution 
for all flux density ranges (in particular, a probability of 
95% is found comparing the intervals at 100-200 mJy and 
> 800 mJy). 

Differences in the 11 distributions for flat-spectrum 
sources can be perceived in Figure |3]for 11 < 1%. However, 
they are probably not intrinsic but induced by instrumen- 
tal effects: at low fluxes (100-400 mJy) only few sources are 
detected with H < 0.1% because the noise rms contribution 
is ~ 0.3%; viceversa, at S > 400 mJy, where the noise is 
practically negligible, a signiflcant tail of very low values of 
n is observed. Finally, we note that artiflcial peaks in the 
distributions at H ~ 0.1-0.25% are produced by the residual 
instrumental polarization. 

The origin of the anti-correlation in steep-spectrum 
sources is not clear yet. A possible explanation is a system- 
atic increase of the mean redshift of sources with decreasing 
the flux density, which would entail a decrease of the Faraday 
rotation measure cx (1 -I- z)~^. A moderate anti-correlation 
between mean redshift and flux density is indeed expected 
based on the T98 model. On the other hand, currently avail- 
able data do not provide any evidence of a positive correla- 
tion between polarization degree and redshift of sources. 

Contrary to M02, we do not flnd any evidence that the 
polarization of flat-spectrum sources depends on the flux 
density. The discrepancy is probably a consequence of the 
different separation in steep- and flat-spectrum sources. In 
M02, the fraction of NVSS flat-spectrum sources is very high 
(~ 44%), and roughly constant at every flux-density inter- 
val. It is possible that, not having taken into account the 
effect of the different angular resolution in the comparison 
between NVSS and GB6 flux densities, M02 have underes- 




Figure 4. The IT distribution for NVSS flat-spectrum sources 
(solid histogram) compared with results of the simulated cata- 
logue (dotted histogram). The dashed curves are the fit given by 
the first term of Eq. 0. 




Figure 5. The same as in Figure|^but for 11 = [0, 1]%. 



timated the spectral index of a significant fraction of steep- 
spectrum sources, misclassifying them as flat-spectrum. 



2.1 The 1.4 GHz H probability distribution for 
flat— spectrum sources 

The Stokes parameters Q and U, measured by the NVSS, 
are the sum of three terms: 
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Table 3. Polarimetric data on extragalactic radio sources. 



Reference 


i/(GHz) 


N. sources 


Comments 


Eichcndorf & Reinhardt (1979^ 


[0.4, 15] 


510 


multi-frequency data 


Tabara & Inoue C 1980 ~) 


[0.4, 10.7] 


1510 


multi-frequency data 


Simard— Normandin et al. fl981) 


[1.6,10.5] 


555 


multi-frequency data 


Perlev (1982') 


1.5,4.9 


404 




Rudnick et al. ('1985) 


[1.4,90] 


20 


fiat; multi— frequency data 


Aller et al. (1992 ~l 


4.8,14.5 


62 


35 fiat, 27 steep; complete sample (Ssghz > 1-3 Jy) 


Okudaira et al. (1993") 


10 


99 


fiat; complete sample {S^qhz > 0.8 Jy) 


Nartallo et al. (1998') 


273 


26 


fiat 


Condon et al. (1998) 


1.4 


~ 2 X 10^ 


complete sample {Si ^qhz > 2.5 mjy) 


Allcr ct al. (1922) 


4.8,14.5 


41 


BLLAC 


^Zukowski ct al. (1999') 


4.75 


154 


28 fiat, 122 steep 


^Lister (2001) 


43 


30 


fiat; 90% complete sample {S^gHz > 1-3 Jy) 


VLA Calibrators 


5,8.5,22,43 


62 


55 fiat, 7 steep 



Q, i7(obs) = Q, U{mt) + Q, [/(noise) + Q, [/(res) , 



(1) 



where Q,[/(int) are the intrinsic values, while the other 
two terms represent the contribution of the noise and of 
the residual instrumental polarization. We find that the dis- 
tribution of Q,U{ohs) for flat-spectrum sources can be re- 
produced adopting the following expression for the intrinsic 
distribution of polarization degrees at 1.4 GHz 1X1,4 (con- 
sidering the flux density from the NVSS catalogue and a 
random distribution of polarization angles): 



with 



/(ni.4 



2.7 + 0.025n5';| 



(2) 



^ ifni.4s;o.i 
o' ifni.4>o.i 



The values of the coefficients (a = 0.51, b = 0.24) are deter- 
mined by the condition that the observed fraction of sources 
with III. 4 < 1% has to be reproduced. The first term in 
Eq. (|5J is the best fit for the Hi. 4 distribution of sources 
with 100 < 5* < 200 mJy and 11 1.4 ^ 1% (it is shown by 
the dashed lines in Figure 2J- The second term corresponds 
to a population of nearly unpolarized sources. In fact, for 
about 20% of the sources, the polarization degree given by 
the N'VSS catalogue is less than 0.4%: these values are only 
upper limits to the intrinsic polarization since they can be 
accounted for by instrumental polarization and noise. For 
such sources we have assumed an intrinsic polarization uni- 
formly distributed between and 0.1%. 

The distributions of the noise and of the residual instru- 
mental polarization for Q and U are assumed to be Gaussian 
with zero mean and variance 0.29 mJy and O.OOIS, respec- 
tively. 

Figures |1| and |3 compare the polarization degree dis- 
tributions of N'VSS sources with the results of our simu- 
lations. The good agreement at every flux-density interval 
is evident and it is confirmed by the Kolmogoroff-Smirnov 
(K-S) test. The probability that the two sets of data come 
from the same distribution is: 6 x 10"^, 0.23, 0.30, 0.60 for 
100-200 mJy, 200-400 mJy, 400-800 mJy, > 800 mJy respec- 
tively. The fit is not satisfactory only for the low polarization 
portion (11 < 0.5%) of the distribution for the 100-200 mJy 
interval (see Figure 0), probably due to our difficulty in cor- 
rectly modelling the instrumental polarization. 



3 POLARIZATION DATA FOR 

EXTRA-GALACTIC RADIO SOURCES AT 

I' > 1.4 GHZ 

Data on the polarization of extragalactic radio sources at 
frequencies higher than 1.4 GHz are rather limited, partic- 
ularly above 5 GHz. At high frequencies, where Faraday 
depolarization should be negligible, data are restricted to 
incomplete samples of few tens of objects (see Table 

Complete samples are essential to avoid selection bi- 
ases. For this reason we have extracted from the MIT- 
Grccn Banlc surveys (Bennett ct al. '1986', Lan gston et alJ 
[1990 ■ Griflith ct al. 1990, Griflith ct al. 1991) a complete 
sub-sample of sources with S^ghz > 1.4 Jy in the follow- 
ing sky areas: 



|6|>10°, -0°.5 < 5 < 19°.5 

5''s;ra<19'\ 20*" < ra < 4*" , 17° < 5 < 39°. 15 

15''.5 < ra < 19'" , 20'' < ra < 2'\5 , 37° < 5 < 51° , 



and we compleme nted it with sources from the Pearson- 
Readhead survey jPearson fc Readheadlfl98ll 1988) at 5 > 
35° and |6| > 10°. Polarization measurements at 4.8 GHz 
were found for 139 (~ 95%) of sources in the final sample 
(hereafter MG-f sample) . Each source was then identified in 
the NVSS and the spectral index between 1.4 and 4.8 GHz 
was estimated. 

Table 2] compares the values of the mean polarization 
degree from samples at different frequencies, while Figures 
|SJ[7|show the polarization degree distributions for fiat- and 
steep-spectrum sources. For flat-spectrum sources, < H > 
increases steadily between 1.4 and 14.5 GHz, although by 
only a factor 1.5, and keeps constant between 15 and 43 GHz 
in the samples of AUor et al. (1993) and of [Lister (200^) 
(note that these flat-spectrum samples are identical, except 
for three sources present at 15 GHz but not at 43 GHz). A 
steady, but weak, increase of the polarization degree with 
the frequency is also found in the case of VLA calibra- 
tors^, simultaneously observed at 5 frequencies between 5 
and 43 GHz. In general, data on flat -spectrum sources indi- 
cate that this population is not strongly affected by Fara- 
day depolarization and that the polarization degree might 

^ http:/ /www. aoc.nrao.edu/~smyers/calibration/ 
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Table 4. Average polarization degree at various frequencies 
for different samples. 



Alli;r t.-L al, 1992 

14.5 GHz 

4.8GHz 



NarUllo el ul, 1996 

273GHz 



5 10 

n 



5 10 15 

n 



Figure 6. fl distribution of flat-spectrum sources for different 
samples at u > 1.4 GHz (solid histograms). Dotted histograms 
show, for comparison, the IT distributions of the same sources at 
1.4 GHz, based on NVSS measurements. 
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Figure 7. Same as in Figure lol but for steep— spectrum sources. 



become frequency independent already at > 15 GHz. On 
the other hand, a significantly higher mean polarization de- 
gree, < n >= 6.6%, corresponding to an increase by nearly 
a factor of 3 from th e value at 1.4 GHz, has been found by 
iNartallo et al.l il998h for a sample of blazars at mm/sub-mm 
wavelengths, where we are probably observing regions very 
close to the active nucleus, whose emission is self-absorbed 
at lower frequencies. Thus, apart from the fact that conclu- 
sions based on incomplete samples need always to be taken 
with care, a large increase of the polarization degree from 
1.4 GHz to 300 GHz can be a consequence of the higher uni- 
formity of the magnetic field in the innermost regions. 

Data on steep-spectrum sources are limited to f < 
15 GHz (see Table I^J, and show a clear increase of the po - 
larization degree with frequency. For the lAUer et al.1 (1199211 
sample < H 1.3%, ~ 2.3%, ~ 4.1% at u = 1.4, 4.8, 
14.5 GHz respectively, while for the MG-I- sample < H > 
changes from 2.4% at 1.4 GHz to 3.3% at 4.8 GHz. The fac- 
tor of 3 increase between 1.4 and 14.5 GHz is comparable 
to the result of M02 from a sample of ~ 130 objects, and 
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CT 


MG-I- 


flat 


85 


1.4 


2.0 


1.9 








4.8 


2.6 


2.1 




steep 


54 


1.4 


2.4 


3.1 








4.8 


3.3 


2.9 


Okudaira et al. ( 1993') 


flat 


89 


1.4 


1.8 


1.7 








10 


2.7 


2.1 


Aller et al. (1992) 


flat 


35 


1.4 


1.8 


1.7 








4.8 


2.3 


2.0 








14.5 


3.0 


2.2 




steep 


27 


1.4 


1.3 


2.1 








4.8 


2.3 


2.4 








14.5 


4.1 


2.2 



Lister (2001) 


flat 


32 


1.4 


2.0 


1.8 






43 


3.0 


2.2 


Nartallo et al. ( 19981 


flat 


26 


1.4 


2.4 


2.0 








273 


6.6 


3.0 



VLA Calibrators 


flat 


55 


1.4 


1.9 


2.1 






48 


5 


2.1 


1.7 






45 


8.5 


2.3 


2.1 






47 


22 


2.4 


2.3 






53 


43 


2.8 


2.4 


Data at 1.4 GHz come 


from the 


NVSS. 









is compatible with a strong Faraday depolarization at few- 
GHz frequencies. 



4 THE n DISTRIBUTION OF 

FLAT-SPECTRUM SOURCES AT HIGH 
FREQUENCIES 

As discussed in Sect.|21 the NVSS data are consistent with 
the polarization degree of flat-spectrum sources being es- 
sentially independent of flux density. If so, the polarization 
properties of bright sources {S > IJy), for which high fre- 
quency data are available, can be assumed to hold also for 
much fainter sources. 

First of all, we check whether the polarization of a 
source at 1.4 GHz is correlated to its value at higher frequen- 
cies. Indeed, the previously mentioned, frequency dependent 
effects (Faraday depolarization and the fact that at different 
frequencies we may eff'ectively observe different emitting re- 
gions) could spoil such correlation. If so, low-frequency data 
would not be representative of high frequency polarization 
properties. The Pearson's linear correlation coefficient r be- 
tween the polariz ation at 1.4 and 4.8 GHz of the sources 
in the MG-I- and lAUer et all ^1992^ samples is ~ 0.6 cor- 
responding to a probability of the null hypothesis, i.e. no- 
correlation, of < 10~^. The correlation is less clear in the 
samples at higher frequencies (r ~ 0.2 and null-hypothesis 
probability of > 10~^). However, these results are probably 
affected by the variability of sources, because we are com- 
paring data taken at different epochs. On the contrary, a 
strong correlation (r > 0.5 and null-hypothesis probability 
< 10~*) between 5 and 43 GHz is found for the VLA cal- 
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Figure 8. Values of e as a function of 11 at 1.4 GHz for MG+ sources (left plot) and for the lOkudaira et al.l ^ 
(right plot). Points with the error bar give < e > with its dispersion. 



sample at 10 GHz 



ibrators, whose polarization is measured simultaneously at 
all frequencies. 

Consequently, the polarization degree of a source at fre- 
quency u can be written in term of its value at 1.4 GHz (Hi. 4) 
as 

n(i/) = ni.4e(ni.4,i^), (3) 

where the factor e is the increase of the polarization degree 
from the value at 1.4 GHz. We find that e itself is a function 
of the polarization degree at 1.4 GHz: Fig.lHl shows e as func- 
tion o f Hi. 4 for the sources in the MG-I- and lOkudaira et alJ 
il993l) samples, while Table |S] gives the average e and its 
dispersion for different bins of Hi. 4. For sources with Hi. 4 
above a few percent, (e) is close to unity at all frequencies, 
but it can take very large values (up to more than 10) if Hi. 4 
is very small. 

We can conclude that: (i) the typical intrinsic polariza- 
tion degree of flat-spectrum sources is of 2-5%; (ii) high 
values of the polarization degree are unlikely (sources with 
H > 10% are rare in the whole range of frequencies anal- 
ysed), probably due to the low degree of the uniformity of 
magnetic fields in radio sources; (Hi) Faraday depolarization 
is probably the cause of the large number of sources observed 
with very low polarization degree (in fact, flat-spectrum 
sources with extreme values of RM h ave b een observed by 
I Stanghellini et al.1 Il998l iTavloil I2OOOI and iPentericci et alJ 
|200(J) . Finally, as already pointed out for < H >, the be- 
haviour of e(14.5GHz) and e(43GHz) is very similar, con- 
firming that the polarization is weakly frequency dependent 
for v > 15 GHz. On the other hand, < e > at mm wave- 
lengths is higher than at GHz frequencies in all Hi. 4 bins. 

Using MG-I- and lOkudaira et al.l il993l) data, we get a 
simple analytical formula for the mean value of e(Hi.4, u): 

< e(ni.4, ly) A{u) exp(-3.2n?:f ) + 0.8 , (4) 

where the frequency dependence appears only in the coeffi- 
cient A{i') which takes the values 38, 50 and 108 at 1/ =4.8, 



10 and 273 GHz respectively. Using these three values, we 
obtain a law, A{iy) = 72 ln(0.75i^°-^ -I- 0.5), which allows us 
to compute < e(Hi.4, u) > at every frequency. The large dis- 
persions around the mean values of e are described by a fit 
similar to Eq. @, obtained interpolating the values of the 
variance of e in the Hi. 4 bins: 

a<,> (H1.4, ly) = B{u) exp(-5H'i':|) -I- 0.5 . (5) 

The factor B is equal to 100 or 210 at u =4.8 or 10 GHz, 
respectively, and it remains constant at higher frequencies 
(we use a linear interpolation at 1/ ^ 10 GHz). 

The above expressions for < e(Hi,4,i^) > and for its 
variance are used to estimate the polarization degree distri- 
butions of fiat -spectrum sources at every frequency in the 
range [5,300] GHz, starting from that of Hi, 4: the increase 
factor e is modelled by a Gaussian distribution which is cut 
at negative e and whose mean value and variance are given 
by Eq. Q and @. In Figure |5| we give examples of the H 
probability distributions between 1.4 and 100 GHz, using a 
simulated sample of 1000 sources [the probability distribu- 
tion at 1.4 GHz is obtained from Eq. ||5J]. The Figure shows 
also the observed distributions ai v ^ 15 GHz, scaled to the 
total number of simulated sources; the good agreement be- 
tween data and simulations is apparent. Moreover, the K-S 
test finds high probabilities (> 0.3) that the observed dis- 
tributions are realizations of the probability distributions 
resulting from our simulations. A significant a greement is 
also fo und with the recent data at 18 GHz by iRicci et alJ 
J2003D . 

Finally, exploiting the T98 model, we generate simula- 
tions of the polarized intensity of the sky at 30 and 100 GHz 
including radio sources and the CMB (see Figure IIOII . We 
can directly compare the importance of the contribution 
of radio sources at the two frequencies: at 30 GHz a lot of 
sources are evident in the map, superposed to the CMB sig- 
nal, while at 100 GHz their presence is far less conspicuous. 
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Table 5. Average values of e and of its variance as a function 
of Hi. 4. 





ni.4 


N 


< e > 


a 


MG+ (4.8 GHz) 


[0,0.1) 


9 


13.5 


13.7 




[0.1,0.4) 


10 


4.52 


4.33 




[0.4,0.7) 


6 


4.13 


1.70 




[0.7,1) 


7 


3.18 


1.84 




[1,1.5) 


11 


1.42 


0.79 




[1.5,2) 


9 


1.88 


0.98 




[2,3) 


10 


0.78 


0.43 




[3,5) 


12 


0.86 


0.45 




> 5 


9 


0.90 


0.46 


Okudaira et al. dggSl flOGHz) 


[0,0.1) 


9 


27.3 


31.0 




[0.1,0.4) 


15 


9.30 


5.38 




[0.4,0.7) 


5 


3.05 


1.45 




[0.7,1) 


8 


2.72 


2.45 




[1,1.5) 


14 


2.72 


2.70 




[1.5,2) 


7 


1.77 


1.48 




[2,3) 


11 


1.30 


0.77 




[3,5) 


14 


0.78 


0.42 




> 5 


6 


0.68 


0.45 



(4.8 GHz) 


[0,0.1) 


1 


1.14 






[0.1,1) 


10 


3.90 


3.28 




[1,2) 


12 


1.77 


1.10 




[2,4) 


8 


0.67 


0.39 




> 4 


4 


1.08 


0.47 


(14.5 GHz) 


[0,0.1) 


1 


44.6 






[0.1,1) 


10 


8.94 


9.38 




[1,2) 


12 


1.98 


0.84 




[2,4) 


8 


1.46 


1.37 




> 4 


4 


0.70 


0.56 


Lister l'2001'> (43 GHz) 


[0,0.1) 











[0.1,1) 


9 


7.43 


4.54 




[1,2) 


11 


2.07 


1.27 




[2,4) 


8 


1.33 


0.93 




> 4 


4 


0.56 


0.41 



Nartallo et al. (1998) 
(273 GHz) 



[0,0.1) 


2 


41.5 


27.5 


[0.1,1) 


6 


8.10 


3.48 


[1,2) 


6 


4.96 


1.95 


[2,4) 


6 


3.05 


1.28 


> 4 


6 


1.10 


0.38 



5 ESTIMATE OF POLARIZATION POWER 
SPECTRA 

Using the statistical characterization of the polarization de- 
gree described in the previous section, we are able to es- 
timate the angular power spectrum (APS) of polarization 
fluctuations due to extragalactic radio sources. We adopt 
the T98 model in order to predict the source number counts 
at cm and mm wavelengths (the model number counts are 
scaled by a factor 0.75 to be in agreement with the WMAP 
data). 

We assume the spatial distribution of sources to be 
Poissonian. The contribution of radio-source clustering to 
the power spectrum is, in fact, small and can be neglected 
if sources are not subtracted down to very low flux limits 
(the clustering is noticeabl e only at relatively sm all fluxes, 
5* < 100 mjy; see T98 and iToffolatti et aljl20o3) . Then the 



Q, 0.2 




0.4 



1.4GHz 

<n>=a.o 




15GHz 

<n>=3.o 




30GHz 
<n>=3.3 




lOGHz 

<n>=2.9 



lOOGHz 
<n>=3.9 



Figure 9. Simulated H distributions at several frequencies be- 
tween 1.4 and 100 GHz. The mean value of H is also indicated. 
The solid line in the plot at 1.4 GHz represents the first term of 
Eq. 0. We report also the observed H distributions, scaled to 
the number of simulated sourc es: 4.8 GHz dat a refer to the MG-I- 
(dashed histogram) and the All er et alj (l99^ (dotted histogram) 
sam ples; 10 GHz and a t 15 GHz data to the lokudaira et al.Nl993l) 
and I Aller et alj ll992l) samples respectively. 



power spectrum is independent of the angular scale, and 
can be computed in the limit of small scales. Under these 
conditions, we have Cei + Cbi = Cqi + Cui iTucci et alJ 
2002). Moreover, because point sources contribute, on aver- 
age, equally to the E-, B-mode and to Q, U power spectra 
jSeliaklll997i. M02), Cee ^ Cse ^ Cgt ~ Cue (in the fol- 
lowing we refer to any of these spectra as the polarization 
APS). 

We estimate the point-source polarization APS in two 
ways: first as the product of the temperature power spec- 
trum, Cit, and the mean squared polarized intensity; sec- 
ond, directly from the number counts of Q (or U). In both 
cases, we exploit the 11 probability distributions obtained in 
the previous section^ 

T egmark fc Efstathioul lll996l) have shown that the APS 
of intensity fluctuations produced by Poisson-distributed 
sources can be written as: 



n{S)S' dS, 



(6) 



where N and n{S) are, respectively, the total and the differ- 
ential number of sources per steradian, and Sc is the min- 
imum flux density of sources that can be individually de- 
tected and removed. 

In analogy to Eq. the Q APS will be 



= N <' 



> 



(7) 



The Stokes parameter Q can be written in terms of the po- 
larization degree p = 11/100 and of the polarization angle in 
the chosen reference system, as Q = 5'pcos(2(j!)). Because 
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Figure 10. Simulated polarized skies at 30 (upper plot) and 
100 GHz (lower plot) generated including extragalactic radio 
sources and CMB radiation (standard CDM model). The sources 
are far less conspicuous at 100 GHz than at 30 GHz. 



we find S, p and to be independent variables in the NVSS 
survey (at least in the case of flat-spectrum sources), we 
obtain 

CQe = N <Q^ > = N < S^p^ cos'(2<^) > 

= N <S^ ><p^ >< cos'(2(^) > 

= 1/2 C« <p^>, (8) 

where the factor 1 /2 arises because of the uniform distribu- 
tion of polarization angles. 

At a frequency u, the mean squared value of the polar- 
ization degree, < >, is: 

nlOO 

v{a,)nldn„ (9) 

'0 

where the probability distribution of II^ at frequency u, 
P(n^), is related to PiUiA) [given in Eq. 

,dni.4 



<pl 10"* 



p(n,) 



p(ni.4)- 



le,(ni.4)-1.12[e,.(ni.4)-0.8]n0;i 



(10) 



Here we have used the relation Il„ = Hi. 46^(111.4) and 
Eq. ^ for the mean value of e^/. For example, we flnd 



< p^ >V^^ 0.041 at 30 GHz and < p^ >V^^ o.046 at 
100 GHz. 

Finally, it is easy to demonstrate that the cross- 
correlation TQ spectrum is 

Ctqi =< S^pcos{2(l>) >= 0. (11) 

The previous method assumes that sources are removed 
from polarization maps using total intensity data. How- 
ever, not all the experiments observing the sky polarization 
are able to measure the total intensity [see, for example, 
the experiments: S POrt dCarretti et aljl2003^ : COMPASS 
jFarese et all2003l) : CAPMAP llBarkatJ20o'^ . In this case, 
sources have to be detected and removed directly from Q and 
U maps. Therefore, we need to compute the Q APS using 
the differential sources counts n{Q), as in Eq. ®: 



Cqi = 



n{Q)Q^ dQ, 



(12) 



where Qc is the lowest value of Q of sources that can be 
individually detected and removed from polarization maps. 
The counts n{Q) can be obtained from the probability dis- 
tribution of the three Stokes parameters that describe the 
linear polarization, 'P{S,Q,U), as 



n{Q) = 



dU / P{S,Q,U)dS. 



(13) 



The probability 7^(5, Q, U) is related to the probability dis- 
tribution of 5", of the polarization degree p = y/ Q'^ + U'^/S 
and of the polarization angle (j) = 1/2 arctan()7/Q), through 



V{S,Q,U) =P{S,p,,p)detJ 



p (j) 
Q U 



(14) 



where J is the Jacobian matrix of the transformation 
{S,Q,U)^ {S,p,<p)- 



detj 



U 



2s^g2 + f/2 



As previously noted, the variables S, p and 
dent, and therefore 



(15) 



are indepen- 



■P{S,P, 



V{S)V{p)V{<j>) = N-^n{S)V{p)/TT; 



(16) 



n(S) is provided by the T98 model and V{p) = lOOP(n) by 
Eq. JISJ. Finally: 



n{Q) = 



dU 



yjWTTP 



s 



dS. 



(17) 



Figure 1111 shows our estimates of the polarization APS us- 
ing the two methods discussed above, at the frequencies 30, 
44, 70 and 100 GHz. In each panel, the two dotted curves 
are computed using Eq. Q (i.e., the first method), and 
Sc = IJy (upper line) or the freq uency dependent de tec- 
tion limit (see caption) obtained bv lVielva et al.l i2003h for 
the Planck mission using the spherical Mexican hat wavelet 
to remove bright sources (lower line). The latter limit, mul- 
tiplied by %/2, is used for Qc (and Uc) in Eq. 1121 (assuming 
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that the source removal algorithm has a similar efhciency 
for polarization maps as for intensity ones). The factor \/2 
arises if the noise in temperature and polarization are re- 
lated by (Jy = (Tp/2. The extrapolation of the polarization 
degree distribution at 1.4 GHz to frequencies between 30 and 
100 GHz is obtained using Eq. Q and the mean value of e, 
Eq. ||1J. The uncertainty on the polarization APS due to 
the high-frequency extrapolation can be estimated by dif- 
ferent realizations of H distributions generated following the 
recipe we describe in section |1] We find that the variance is 
small, a few percent of the average spectrum (the number 
of sources is large enough to assure the sample variance to 

be negligible). 

A s discussed in section |1] if the data of iNartallo et all 
l)l998h are disregarded, the polarization degree appears to 
be weakly dependent on the frequency for v > GHz. As- 
suming the H distribution to be constant above 15 GHz, the 
amplitude of the polarization APS decreases weakly with 
respect to previous estimates: ~ 10% at 30 GHz and ~ 30% 
at 100 GHz (the differences are visibly appreciable only in 
the 100 GHz panel of Fig. lllH . These latter estimates can be 
considered as lower limits and, together with the ones us- 
ing the high-frequency extrapolation, provide us with two 
boundaries within which we expect to find the real polariza- 
tion APS. 

We have verified that the results from the two methods 
agree when the bright source removal is equivalent. Other- 
wise, as illustrated by Fig. 1111 estimates based on Eq. ^ 
are significantly lower than those using Eq. 1121 : the differ- 
ence is nearly a factor 10 if we take Qc ~ Uc ~ V^Sc (in 
the real case, however, we expect that Qc = Uc > V^Sc, 
further increasing the difference) . This is easily understood: 
adopting the limit on total flux density we are subtracting 
many more sources than we can do using only polarization 
maps. 

The previous estimates are obtained assuming that the 
polarization degree distribution for steep-spectrum sources 
is the same as for the flat-spectrum ones, although the re- 
sults reported in Table |3] and by M02 indicate that steep- 
spectrum sources are, on average, more polarized, at least 
at low frequencies. However, the contribution of steep- 
spectrum sources to the polarization APS is anyway very 
small at i/ > 40 GHz. 

The likely correction to our results at 30 GHz can 
be estimated as fol lows. The VSA survey at 15 GHz 
ijWaldram et al.ll2003h flnds a 25 percent of steep-spectrum 
sources for S15 ^ 100 mjy and a 44 percent for S15 ^ 
25mJy, corresponding to a contribution to Cie of 30% if 
Sc ~ IJy and of 43% if Sc = 200 mJy (we assume that 
the fraction of steep-spectrum sources increases to 75% if 
5*15 < 25mJy). Their contribution to the APS is reduced 
by a factor 3 at 30 GHz if their average spectral index is 
—0.8 between 15 and 30 GHz (a rather conservative as- 
sumption since the high frequency spectra steepen due to 
electron ageing effects). From the NVSS catalogue we find 
that < Hf 4 >^^^~ 2.6% for steep-spectrum sources with 
5i.4 ^ 200 mJy. Now, assuming that the polarization degree 
increases, on average, by a factor of 3 from 1.4 GHz to high 
frequencies, we obtain that our estimates of the polarization 
APS by Eq. must be multiplied at 30 GHz by 1.3 or 1.4 in 
the case that S'c = IJy or 200 mJy respectively (dot-dashed 
lines in the 30 GHz plot of Figure lllll . 



6 DISCUSSIONS AND CONCLUSIONS 

The main goal of the present paper is to provide estimates 
of the contamination of CMB polarization maps by extra- 
galactic radio sources. This task is particularly complicated 
because of the lack of polarimetric data at the frequencies 
where the CMB is observed. 

The richest data on the polarization properties of ra- 
dio sources comes from the NVSS survey, which provides a 
very large and complete catalogue of extragalactic sources at 
1.4 GHz. Exploiting the spectral information obtained com- 
bining the NVSS with the GB6 survey at 4.8 GHz, we have 
analyzed the properties of steep- and flat -spectrum sources. 
The former population shows an anti-correlation of the po- 
larization degree with the flux density, already pointed out 
by M02. On the contrary, for flat -spectrum sources the dis- 
tributions of the polarization degree for different flux-density 
ranges do not show any signiflcant variation and its mean 
value is constant and quite low, around 2%. For this class of 
sources we flnd a flt that accurately describes the polariza- 
tion degree distribution for all sources with 5* ^ 100 mJy. 

Exploiting the available data ai 1/ > 1.4 GHz, we have 
investigated how the polarization degree of sources varies 
with frequency. Data on steep-spectrum sources, available 
only up to 15 GHz, highlight a strong increase of the mean 
polarization degree with frequency, consistent with substan- 
tial Faraday rotation measures. On the other hand, the mean 
polarization degree of flat-spectrum sources shows only a 
weak increase between 1.4 and 15 GHz, and appears to re- 
main essentially constant at higher frequencies. However, the 
high values of the polarization degree found at mm wave- 
lengths by Nartallo ct al. ( 1998) may be an indication that 
new, more polarized, components show up there. 

For flat-spectrum sources, we have derived an analytic 
function allowing us to extrapolate to any frequency the dis- 
tribution of polarization degrees at 1.4 GHz. We notice that 
only sources with Hi. 4 < 1% increase signiflcantly their po- 
larization degree between 5 and 43 GHz. These sources may 
be the only ones with really strong Faraday depolarization at 
1.4 GHz. In general, at v > 5 GHz the polarization degree of 
flat-spectrum sources is typically of a few percent (it rarely 
exceeds 10%), and is weakly dependent on the frequency. 

Using this analytic function, we estimate the power 
spectrum of polarization fluctuations induced by extragalac- 
tic radio sources, by means of the two methods described in 
Section |S] In particular, we consider for the flrst time the 
case whereby the removal of the brightest sources can be 
only done using polarization data. In this case, the source 
subtraction is far less efficient than if total fiux data can 
be used and, as a consequence, the amplitude of the source 
power spectrum can be up to a factor of 10 higher. Total in- 
tensity data are, therefore, essential information in order to 
reduce the contamination of CMB polarization maps. In Fig- 
ure llll we show the estimated contribution of undetected ra- 
dio sources to the polarization signal observed by the DASI 
experiment at 30 GHz: if indeed all the sources with the fiux 
density higher than 50mJy have been removed, the residual 
contamination is extremely small. 

We have focused our analysis on flat-spectrum sources, 
since they are the dominant population at least up to 
100 GHz. Steep-spectrum sources are taken into account 
only at 30 GHz, while they are disregarded at higher fre- 
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Figure 11. Polarization power spectrum produced by extragalactic radio sources compared to the E— and B-mode CMB spectrum (CDM 
model with Qa = 0.7, re-ionization optical depth r = 0.17 and tensor-to— scalar ratio r ~ 0.5). The short-long dashed lines are the 
j?— m ode component induced by gravitational lensing. The CMB and lensing spectra are computed with CMBFAST iSeliak &: Zaldarriagal 
Il99^. The APS for extragalactic radio sources are computed using the first method with Sc = 1 Jy (upper dotted lines) and Sc estimated 
bv lVielva et al.l <2003l) for the Planck mission (200, 400, 250, 200 mjy at 30, 44, 70, 100 GHz; lower dotted Unes), and using the second 
method with Qc = Uc = V^Sc, where Sc is that from Viclva ct al. (2003} (dashed lines). In the plot at 30 GHz, the dot-dashed lines 
corresponds to the estimates by the first method taking into account the correction for the contribution of steep-spectrum sources (see 
text); the long dashed line is our estimate for the contribution of undetected sources in the DASI experiment. In the plot at 100 GHz, 
the lower dashed line is the APS computed assuming the H distribution to become frequency-independent for u > 15 GHz. 



quencies. About the GPS sources, large uncertainties are 
still present on their abundance. The WMAP data at 33 GHz 
detected 16 sources with inverted spectra (q < —0.4) over 
a total of 155 source s with S > 1.2 Jy. The multifrequency 
data summarized bv iTrushk in (2003) indicate that most of 
them are likely to be blazars (flat-spectrum quasars or BL 
Lacs) whose radio emission is dominated by a single emitting 
region (a knot in the jet) caught during a flare. Thus, the sur- 
face density of bright GPS sources peaking at high radio fre- 
que ncies is likely to be se veral times lower than predicted by 
the lDe Zotti et alJ i200(f ) models, even in the case of an in- 
trinsic distribution of peak frequencies not extending above 



~ 200 GHz. According to the latter model (which, in the 
light of WMAP sources, provides rather generous upper lim- 
its), at u — 100 GHz, the contribution of GPS sources to the 
total intensity APS is at the few percent level. Furthermore, 
the polarization degree of bona-fide GPS sources is very low 
at cm wavelengths. GPS galaxies generally have fractional 
polarization below 0.3% (many are undetected below 0.1% 
levels); GPS quasars have higher polarization t han galaxies, 
but lower (at 6 cm) th an non-GPS quasars l! 0'Dea et al] 
ll990 tlStanghellini et aD fl998: StanEhclhni 1999). Thus, the 
contribution of GPS sources to the polarization APS is likely 
to be even lower than to the total intensity APS. 
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It is interesting to compare our results on extragalactic 
radio sources with the E- and B-mode spectra of the CMB 
radiation (see Figure 111(1 . First of all, we notice that the 
CMB E-mode APS increases with ^ in a similar way to the 
point-source spectrum in the range 30 < £ < 10^. At these 
angular scales, radio sources yield a relevant contribution to 
the E-mode polarization only for < 40 GHz; in any case, 
also at these frequencies, they are not a serious problem for 
measurements of the CMB E-mode APS when the brightest 
sources are removed. At higher frequencies the cosmological 
signal is dominant at least up to £ ~ 2000 or more. 

On the other hand, extragalactic radio sources can be 
a critical factor for the detection of the CMB B-mode com- 
ponent: the cosmological Cbi is characterized by a peak at 
I ~ 100 whose amplitude is directly related to the inflation 
model and, in particular, to the ratio between the ampli- 
tude of tensor and scalar perturbations. It is well known 
that an important limitation to the detectability of such 
peak comes from B- modes induced by gravitational lensing 
llKnox fc Song l2002h . especially on scales I > 100. How- 
ever, our analysis highlights that, at least for f < 100 GHz, 
B-modes produced by extragalactic radio sources are even 
more critical. A better understanding of polarization prop- 
erties of dust in the Milky Way and in external galaxies 
is necessary to establish whether frequencies higher than 
100 GHz may be more suitable to investigate the B-mode 
polarization. 
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